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Abstract 



■ A gravity theory called scalar-tensor-vector gravity (STVG) has been recently developed and 



succeeded in solar system, astrophysical and cosmological scales without dark matter [J. W. Moffat, 
J. Cosmol. Astropart. Phys. 03, 004 (2006)]. However, two assumptions have been used: (i) B{r) = 



A~^{r), where B{r) and A{r) are goo and in the Schwarzschild coordinates (static and spherically 
symmetric); (ii) scalar field G = Const, in the solar system. These two assumptions actually 
imply that the standard parametrized post-Newtonian parameter 7 = 1. In this paper, we relax 

^ : .™ ... S.VO ... . _ <P. 

approach. With abandoning the assumptions, we find 7 7^ 1 in general cases of STVG. Then, a 

o . 

On . version of modified STVG (MSTVG) is proposed through introducing a coupling function of scalar 



field G: 0{G). We have derived the metric and equations of motion (EOM) in IPN for general 



^ ■ matter without specific equation of state and N point masses firstly. Subsequently, the secular 

periastron precession tu of binary pulsars in harmonic coordinates is given. After discussing two 
PPN parameters (7 and (3) and two Yukawa parameters (a and A), we use cu of four binary pulsars 
data (PSR B1913+16, PSR B1534+12, PSR J0737-3039 and PSR B2127+11C) to constrain the 
Yukawa parameters for MSTVG: A = (3.97 ± 0.01) x lO^m and a = (2.40 ± 0.02) x 10"^ if we fix 
|27-/?-l| =0. 
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I. INTRODUCTION 



With tremendous advance in the accuracy of observations, Einstein's general relativity 
(GR) has passed nearly all the tests in the solar system. However, alternative gravity 
theories still stand up for explaining some exotic phenomena such as the dark matter and 
for possible violation of general relativity in the future higher precision experiments. Among 
them, Moffat p.] proposed a scalar-tensor-vector gravity (STVG) theory in which there are 
three scalar fields and a vector field besides the metric tensor. Three scalar fields are 
respectively the scalar field G that origins from the Newtonian gravitational constant, the 
coupling function of the vector field u and the rest mass of the vector field m which controls 
the coupling range. They all change with space and time. The vector field 0^, which is 
associated with a fifth force charge, corresponds to the exchange of a massive spin 1 boson 
and couples to the ordinary matter. Through introduced 0^, a Yukawa-like force was added 
to the Newtonian inverse square law. This leads to a satisfied fit to galaxy rotation curves 



and the Tully-Fisher law 
cosmological observations 



2|. Besides, the theory has been_^used successfully to explain 
3|, the motion of galaxy clusters 



, the Bullet, Cluster Q, the 



velocity dispersion profiles of satellite galaxies [10|| and globular clusters {ll] without exotic 
dark matter. By studying STVG attentively, we found it uses two assumptions in STVG: 

1. The metric components goo and Qrr in a spherically symmetric field have the following 
relationship: B{r) = A^'^{r) , where B{r) and A{r) are goo and Qrr respectively in a 
Schwarzschild-like coordinate system; 

2. For the solar system, the running of gravitational constant G is zero, namely, G = 
Const. 

With the above assumptions, let us analysis what results they lead to. A general form of 
the metric in the Schwarzschild coordinates reads 

ds^ = -Biry^dt" + A(r)dr2 + r'^d?^, (1) 

where d^VL = d6'^ + sin^ Odip"^. Comparing Eq. ([T]) with the following parametrized post- 
Newtonian metric, 

ds' = - (l - c'df + 27^) dr' + rM'Q, (2) 



within the post-Newtonian precision C(l/c^), it is clear that 7 = 1 under assumption (i). 
On the other hand, since only the scalar field G affects 7 in STVG (for details see subsection 
D of section II), assumption (ii) gives 7 = 1 directly. In summary, these assumption impose 
a constraint on the theory: the standard parametrized post-Newtonian (PPN) parameter 7 
equals to 1. Therefore, these two assumptions are reasonable and consistent with the current 
measurement of 7. 

However, we prefer to a general form in natural science without any imposed limitations. 
For the first assumption, we have no strong reason to take it. For the second assumption, 
the search for the Newtonian gravitational constant never stops. For instance, planetary 
and spacecraft ranging, neutron star binary observations, paleontological and primordial 



12|. This is 



nucleosynthesis data allow one to constrain the variation of G with time 

the reason why we relax these two assumptions and study STVG further. Firstly, Let 

us investigate the numerical value of 7 in STVG attentively with abandoning these two 
assumptions in comparison with Brans-Dicke theory. We expand metric as follows 

^700 = -H-e2iV + e^L + 0(5), (3) 

go^ = e%, + 0(5), (4) 

Qij = Sij + t^H.j + 0(4), (5) 



by using Chandrasekhar's approach [13|], we have found and Hij only depend on the scalar 
field G and matter in STVG (for details see subsection D of section II) and where e = 1/c 
and 0{n) means of order e". Then, we can compare STVG with Brans-Dicke theory (BD) 
Q. The action of BD is 

SBD = ^j {ci>R - ^09'"^^ V^d'x + Sm, (6) 
where <^o is a coupling constant {uq usually is used in BE*, here we use <^o to avoid confusion 



with another symbol in the context). For STVG 
the metric and G is 



15l | , the action by only considering 



Comparing Eq. ([6]) with Eq. (JTj), we have G = 1/0 and ^0 = —1/2. In BD, 7 = (^o + 
l)/(<Jo + 2) and when <jo — > C)0, 7 = 1. With the correspondance between BD and STVG, we 
obtain 7 = (^0 + l)/(^o + 2) = 1/3 7^ 1 in Moffat's STVG. (But in some Refs.Q, S 3, the 



value of ^0 is 1/2 due to the sign of the scalar field action changed and correspondent 7 is 
3/5.) In subsection D of section II and Appendix A, we give a strict proof about 7 7^ 1 of 
STVG in general cases. 

Although STVG is a theory that violates EEP, with relaxing the above two assumptions, 
the departure of 7 from 1 should not be out of the range restricted by current experiments. 
7or example, the measurement of 7 in the Cassini experiment gives 7 — 1 = (2.1 ±2.3) x 10~^ 
1^. The data analysis for this result is conducted under standard PPN framework. One 
might argue that STVG modifies the Newtonian law by introducing a Yukawa potential, 
which breaks the PPN framework and therefore the constraint given by Cassini may not 
be used directly here. Another measurement of 7 comes from the Lense-Thirring effect 
(precession of a gyroscope located near a rotating massive body). Lense-Thirring effect 
will cause advance of the ascending node of the orbit of an earth satellite, which depends 
only on a force perpendicular to its orbital plane. But the Yukawa force is confined in the 
orbital plane so that the measurements of Lense-Thirring effect can give a clear and direct 
constraint on 7 in STVG. Although the Lense-Thirring effect for LAGEOS and LAGEOSII 
are at a precision level of 10% 17|], worse than the Cassini experiment, this experiment 
demands that 7 can not departs from 1 too much (We will explain this in subsection A of 
section III in detail). This is the main motivation to propose a modified STVG theory in 
this paper when we abandon these two assumptions. It is more competitive for present and 
future experiments. 

With discovery of the binary pulsar PSR B1913+16 by Hulse and Taylor [l^, binary 
pulsars promises an unprecedented opportunity to measure the effects of relativistic gravi- 
tation (see Q, [2^ for a review). For example, pulsar timing has provided indirect evidence 
for the existence of gravitational wave s |2ll| . the binary pulsars data can constrain the ex- 



istence of massive black hole binaries [22|], and the binary pulsars can a. 



so test the effects 



23|. In addition, bi- 



of strong relativistic internal gravitational fields on orbital dynamics 
nary pulsars could help us to test various gravity theories. By fitting the arrival time of 
pulsars, observational parameters of binary pulsar are obtained in high precision. They are 



"physical" parameters, "Keplerian" parameters and "post- Kepler ian" parameters 2J, |25| 



A very important class is "post-Keplerian" parameters |26l . l27l | which include the average 
rate of periastron shift u, redshift dilation 7, orbital period derivative Pb, and two Shapiro- 
delay parameters s and r. It worthy of noted that the periastron advance for binary pulsars 
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could reach several degrees per year, which is about 10^ more than the perihelion advance of 
Mercury. Hence, the relativistic effects from binary pulsar are more remarkable than other 
celestial systems. We will take a sample of binary pulsars for testing MSTVG by using 
their uj. In this paper we chose four best studied pulsar binaries. They are respectively 
PSR B1913+16, PSR B1534+12, PSR J0737-3039, and PSR B2127+11C, which are double 
neutron star binaries. We mainly focus on io of these four binaries data to constrain the 
parameters of MSTVG. 

Through introducing a coupling function of the scalar field G: 0{G), we obtain the metric 
and equations of motion (EOM) for general matter without specific equation of state and 

point masses in the first order post-Newtonian (IPN) approximation by Chandrasekhar's 
approach. Then, the secular periastron precession of binary pulsars for IPN in harmonic 
coordinates is derived. After two PPN parameters (7 and (3) and two Yukawa parameters 
(a and A) are discussed and compared, we fix 7 = 1 and P = I and constrain the Yukawa 
parameters. This constraint coming from binary pulsars systems on a and A are consistent 
with the results from the solar system such as the earth-satelite measurement of earth gravity, 
the lunar orbiter measurement of lunar gravity, and lunar laser ranging measurement to 
constrain the fifth force. 



In what follows, our conventions and notations generally follow those of [28|, the metric 
signature is (-, +, +, +). Greek indices take the values from to 3, while Latin indices 
take the values from 1 to 3. A comma denotes a partial derivative, and semicolon denotes a 
covariant derivative. Bold letters denote spatial vectors. The plot of this paper is as follows. 
In the next section, the equations of motion and u for binary of MSTVG in IPN are given. 
Subsequently, in the third section, we discuss parameters of MSTVG in detail. We then 
constrain the parameters of MSTVG by means of fitting u for four binary pulsars and deal 
with the details of discussion about results in the forth section. Finally, constraints method 
and results are outlined in the last section. 
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II. THE THEORY 



A. Action and field equations 



Through introducing a couphng function of the scalar field G: 6{G), we adopt the fol- 
lowing modified action for STVG 

(8) 



where denotes all the matter fields. In Eq. ([8]), the matter fields \E' interact with both the 
metric field and the vector field. This means that the trajectory of a free-fall test particle 
depends not only on the space-time geometry but also on the vector field so that it violates 
the Einstein equivalence principle (EEP). Although the current experiments verify EEP to a 
very high accuracy in the Solar System scale, violations of EEP at galactic and cosmological 
distance scales can not be ruled out. In Eq. ([8]), the Lagrangian densities of the gravitational 
field, vector field and scalar fields are respectively 



C^ 



16nG 
1 

UJ 

c 



(i? + 2A), 



4 



^m'e<p,r - n(0) 



c 

IQttG 



+ UJ-^UJ-u 1 - 



(9) 
(10) 

(,11) 



where A denotes cosmological constant and Vx{X) denotes the self- interaction potential 
function of an field. Besides, k = /h? and B^j^^ = 4>u,fi — (pti,u- Where c is the ultimate speed 
of the special theory of relativity and h is the reduced Planck constant. It is noted that in 



Moffat's STVG 



15|, OiG) 



-1 in Eq. ([TT]). 



When we omit A and all the self- interaction potentials Vx{X), equations of gravitational 
field are obtained by variation of the action ([8]) with respect to g'^^. 



R 



8nG 



T, 



V + -^'^ I B^j_f^BJ^ — -g^uB^xB 



G 



4 



1 fG^ 



r) G-i^G'^ 



where Tnv is the stress-energy-momentum tensor of matter which is defined by 



5S 



M 



5g 



(12) 



(13) 



Following [291] and [30|, l3l|, we define mass, current, and stress density as 





(14) 




(15) 


a,, = c^T^. 


(16) 



It is worth emphasizing that, in these definitions, the matter is described by the energy- 
momentum tensor without specific equation of state. 
Variation of the action with respect to (p^ yields 



where is a "fifth force" matter current defined as 

5Sm 



'-9 — 
c 



(17) 



(18) 



We further define 



J"" = J°, T = eJ\ (19) 

When we substitute k = c^/h? and B^,^ = — (p^^u into Eq. (fT7|) . we can obtain the Proca 
equation if u; = Const, and J'^ = 0. 

Variation with respect to the scalar fields yield respectively 

(9(G) H- 3)iG-. - 2(9(G) + 3)^ - 



8nG 



1 

—m'.„ - 
m 



T -uo- 



:m' m-i, — — — G' m-y 



Gm 



:G'''UJ._„ 



''^ G' 

B. Perturbation of MSTVG 



(20) 
(21) 
(22) 



Following the approach in 13|,|32|, we deal with MSTVG in the form of a Taylor expansion 
in the parameter e = 1/c, similar to the expansion of the metric tensor in Eqs. ([3]), (jl]) and 
([5]). For the expansions of vector field, 



(0) 2(2) 

(1) 
eLfi. 



(23) 
(24) 



and for the expansion of scalar field G 



G = 6-0(1 + = G'ofl + e'G + e^G J, (25) 

d{G) = eo + e^^ + h2e + ■■■ , (26) 

(2) /I (2)2 (4)\ 

e{G) = Oo + e^OiG + e^i-e2G + OiG l (27) 

^^^(^^^^^^^^ 

where Go is the background value of scalar field G. The expansions of others scalar fields 
are 

m = mQ ^ + e'^m + e^rrtj , (29) 
Lu = iUo(^l + e^S + e^S^ , (30) 
where ttiq and ujq are respectively the background values for the scalar fields m and uj. 

C. Gauge condition 

We u.e the gauge »ndition imposed on the component of the metric tensor proposed by 
Kopeikin & Vlasov ^] as follows: 

Go 



^.^9''] =0. (31) 



33|. 



Noted that the scalar field G in MSTVG is in the inverse ratio to the scalar field in 
To IPN order, this gauge gives 

o/l 1 (2) \ 

\2^'' ~ 2^'' ~ ~ ^'V " °' ^^^^ 

and 

/ 1 1 (2) \ 

lLk,k - ^H,t - -Nt + gA = 0. (33) 
Through covariant divergence of Eq. f|T7|) , we derive a useful formula 
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D. First order post-Newtonian approximation for MSTVG 



Based on fields equations of Eqs. ([T2]), ([T7D, ([20D, ([2lD and ([22]), we obtain tlie result of 
MSTVG in IPN by using the gauge conditions ([32]), (I33D and Eq. ^ as follows. 

(2) 

The equation for N and G are 

UN = -Svr^cr, (35) 

(2) 

□ G = (1 - 7)47r^a, (36) 

where □ is the D'Alembert operator in the Minkowski spacetime and Newton's constant Q 
is related to the constant Gq by 



Metric Hij is 



DHij = -8-f7iga5ij, (38) 



where we can define Hij = 6ijV and 

Appendix A gives mathematical details of the derivation of this important formula. A 
rigorous presentation in Appendix B identifies the parameter 7 defined in Eq. ( I39l) is just 
the PPN parameter with the same symbol. 

With 9o = —1, we obtain 7 = 1/3 which just reduces to STVG. Then, Eq. fl37j) becomes 

2 



. ^^-G„. (40) 

It is evident that Q = Gq when 7 = 1. The above shows that only the scalar field G enters 
the metric and Hij and we have given a strict proof about 7 7^ 1 in general cases of STVG 
by Chandrasekhar's approach, (see Appendix A and B) 
Other metric components in IPN are respectively 

DLi = 8(1 + 7)vre?ai, (41) 



□L = -AttG 



(37 - 2/3 - l)Na - 2(1 - 7)(Tfcfc 



2 



+ (37 + l)^JV° + ^(l + 7)^oA^^'^o 

(0) . (0) " 



^(3AN', (42) 
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TABLE I: PPN parameters in MSTVG. 



7 


/9 




Oil 


0-2 


«3 


Ci 


C2 


Cs 


C4 


eo+2 


^ + 2(eo+3)(eo+4)2 



























where A = and 



= l + . (43) 

^ 8(1+7) 2(«„ + 3)(9„ + 4)2 



Appendix B identifies (3 in Eq. (H3l) as the corresponding PPN parameter with the same 
symbol. In that Appendix all the other PPN parameters except 7 and (3 are proved all zero 
for MSVTG. The ten PPN parameters for MSVTG are listed in Table E 
Equations for the vector and other scalar fields are 

^(0) 9,9(0) J° ,,,, 
□Vo-mofcVo = , (44) 

dm = 0, (45) 

(2) 

□ J = 0. (46) 

In Eq. fl44p . if we consider a vacuum case (J'^ = 0), we find that the speed of vector field 
is not equal to the velocity of light when we assume a plane wave solution which may cause 
chromatic dispersion in vacuum. 



E. Equations of motion for MSTVG in IPN 



Based on Ref. [3J], the equations of motion (EOM) derived from T.^'^ = is equivalent 
to Ci^j^ = if EEP is satisfied. However, STVG violates EEP due to the vector field. So, 
equations of motion in MSTVG have to be derived from the covariant divergence of the 



Einstein tensor G^^*^ = 0. On the other hand, according to 



34( 1 . we can infer that the 



equations of motion in MSTVG must include the contribution of the vector field besides the 
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matter and metric. The Einstein tensor G^'^ in MSTVG is 



2^ 



SttG 



UJ 



Xk 



4^ 



m^k^ 1 



2^ 2 



Then, the Bianchi identities G!^ = yields the momentum equation 

^ AT , 70 



+6^ -aVN. + aL:.t - -fxL , 



5 1 

5 111 
+ 2^jfcKfc ~ + ^CTkkN^i - -akkVi 

T/^,.) 7O (1) TO, (2) .0 

Tk W Tk « , TO , Tk 



+Vk,J' - V.,kJ'' - VAJ't + J%) 1 = 0. 

G"^ = yields the continuity equation 



(47) 



(48) 



(49) 



F. N-body pointlike for MSTVG in IPN 



Considering an N- 
notation adopted by 



3ody system of nonspinning point masses for simplicity, we follow the 



35l | and use the matter stress energy tensor as follows 
c2T^'^(x,t) = 5^/i,(t)t;X5(x - Yait)), 



(50) 
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where 5 denotes the three-dimensional Dirac distribution, the trajectory of the ath mass 
is represented by ya(t), the coordinate velocity of the ath body are = dyait)/dt and 
= (c, Vq), and fj,a denotes an effective time-dependent mass of the ath body defined by 

f^ait) = ( , (51) 



where subscript a denotes evaluation at the ath body and Ma is the constant Schwarzschild 
mass. Another useful notation is 

,2" 



1 + ^ 



(52) 



where = v^. Both fia and fla reduce to the Schwarzschild mass at Newtonian order: 
fia = Ma + 0{e^) and Jia = Ma + 0{e^). Then the mass, current, and stress densities in Eqs. 
([11]), f|T5|) and f|T6l) for the N point masses read 



= ^/ia5(x- ya(t)), (53) 

a 

= 5^^X5(x-y,(t)), (54) 

a 

(Tij = ^^a^^X5(x-y„(t)). (55) 

a 

We now assume a "fifth force" matter current 

Cj''{x,t) = J2QaV^6i^ - ya(t)), (56) 



where the "fifth force charge" is Qa = /«v Gujjj,a{t) which coupled with ordinary mass through 
coupling function uj, where k is a dimensionless constant. This assumption implies that this 
charge is proportional to ordinary mass as in [l^. Substituting Eq. (!56l) into Eq. (fT9|) . we 
obtain 

jO = J]«:yG^/i,5(x-y,(t)), 

a 

f = 5^«:VGcu/i,t;^5(x-y„(t)). (57) 

a 

Then, we obtain the following form in Newtonian approximation, 

Ar = 2A-i{-47re;a} = 2^^ + 0(2), (58) 



a 



13 



and 



(0) 
^0 



1 



K 

An 



-mofclx— z 



X — Z 



-J°(z,t)d=^z 



:i + 7)g ^ M,e— o^-" ^ ^^^^^ 



(59) 



Through integration of Eq. (148!) . we have EOM in Newtonian approximation as follows 



dt 



E ^---i 1 - !^^^f±^e-™'- ( 1 + ,n„kr^ \ + 0(2). (60) 



From Eq. (160|) . the gravitational potential in Newtonian approximation for MSTVG is that 



(61) 



On the other hand, Fischbach and Talmadge [36] provided the following potential 



ae ^ 



(62) 



which includes a usual Newtonian gravitational potential and a Yukawa-type one. Compared 



between Eq. (161!) and Eq. (1621) . we define uiok = 1/A and a 
becomes 



K^(l+7) 



cjQ. Then, Eq. ([60]) 



dt>! 



(63) 



If MSTVG could explain galaxy rotation curves without exotic dark matter, it must have 
a > from Eq. ( l63l) . 

The next step is to work out and V easily in IPN approximation as 



'a I 'a 

a ^ a 



a b^a 



Avl - {UaVaf 

1 + ae"'""/^ 



+ 2(2-37)5^5^ 



1 + 



rab 



0{3), 



G^MgM, 

TaTab 



V ^ 2,J:^ + 0(2), 



by the relation of fia that 



(64) 
(65) 



b^a 



+ 0(4) 



(66) 
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where = |x - y^l and rat = \ya - Ybl- For Li 



i. = -2(l+7)E^< 



(67) 



Because a in the Newtonian order is very tiny (see Table UTTl) . we omit any couphng terms 
in the magnitude of ae^. Thus, it yields 



a b 



ran 



At last, by integration of Eq. (gH]), we obtain EOM for MSTVG in IPN: 



dt 



~"'ab ^ / , „2 "-ab 



b^a 



ab 



b^a 
?2 i\j2 



+6^ 2(7 + ^-^<b + (27 + 2/5 + 1) 5^ : ' n 



bjta 



bj^a 



b^a 



7t;^ - (1 + 7)t;2 + -{nabV^f + 2(1 + 7)(t;at;b) 



n 



ab 



2(1 + 7)(nabt;a) - (27 + l){nabVh) 



ab 



bf^a c^a,b 
1 

2 



ab 



1 

Tbc 



{2(3 - 1)— + 2(/3 + 7)— - ^(n,,n,,) 



fabric 



(68) 



(69) 



b^a cj^a,b 

where n^j, = (y* — yl)/rab, ^ab = Va — v^, and scalar products are denoted with parentheses, 

1 return to Einstein-Infeld-Hoffmann (EIH) equations 



e. g., [nabVab) = T^ab-^ab- Eq. (|69 
of motion in the PPN formalism 
field in MSTVG. 



wi 



24 



371] when we eliminate the contribution of the vector 



G. Secular periastron precession of binary pulsar for MSTVG in IPN 



When we consider only two body (Mi and M2) in Eq. ( 169|) . EOM of a point-mass binary 
for MSTVG in IPN yields 

r = a^v + aiPTv, (70) 
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TABLE II: Summary of the PPN parameters 7 and /? in several gravity theoties. 



Parameter GR 



ST 



MST 



STVG 



MSTVG 



M TeVeS 



7 



UJQ + 2 
1 I m 

-'"^(2a;o+3)(2cJo+4)^ 



1 - 



1 + 



2a2 



1 + 



60+2 

eo+4 
61 



2(eo+3)(9o+4)2 



Violation of No 
EEP? 



No 



Yes 



Yes 



Yes 



No No 



where 



^IPN 



GM 



r \ r 
l + a\ 1 + - jexp(--) 



n, 



QM r 



1 



7 + 31/ t,'^- 2 7 + , 9 + V 



I r 



3 .2 

—ur 

2 



(71) 

n-2( 1 + 7 - z/ )rv [>(72) 



where M = Mi + M2, r = (^12^12), v = vi — V2, v = MiM2/M^, r = ri2, n = ni2 and 



a = dvi2/dt. Compared with [38!], additional term comes from the Yukawa force, 



namely. 



the contribution of the vector field in MSTVG. By the aid of the averaging method 391], the 
secular periastron advance for a binary pulsar in IPN is 

du 1 nap^ 



-/^ + (2 + 27-/3)^^ 



(73) 



dt 2 A2 p 

where n^o? = QM, p = a{l — e^), a is the semi-major axis and e is the eccentricity of the 
binary. Here the periastron shift of Yukawa force is 

du 1 nap"^ 



-p/X 



dt 



Yukawa 



2 A2 



(74) 



which is different from the result given by 401]. In this paper, we keep the angular momentum 
in the process of using the averaging method instead of approximate treatments such as 
exp(— r/A) ^ 1 — r/A |40] to calculate the u. 



III. SUMMARY OF PARAMETERS FOR MSTVG 

There are four parameters in MSTVG. They are two PPN parameters (7, (3) and two 
Yukawa parameters (a. A). In this section, we mainly discuss these four parameters. 
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A. Review and discussion on PPN parameters 



For a slow motion and weak field, PPN formalism introduces 10 parameters in a post 



Newtonian metric to include various gravity theories 
back to Eddington-Robertson-Schiff formalism 



41 



43 



44 



42| . PPN formalism can be traced 



45l ] , which introduce two PPN param- 



eters 7 and (3 by only considering one point mass. At the same time, these two parameters 
were endowed with some kind of meaning. For example, 7 denotes the level of space cur- 
vature and (3 can be treated as described the nonlinearity in the superposition law 24 1. 



Tab le [Til lists the results of 7 and (3 in: (1) General relativity (GR); (2) Scalar- Tensor theory 



46|, CkST, where luq and Ui are coefficients of the coupling function of scalar field 



e(0) = ^o+^i^+0(a if- 



1 + 0); (3) Multi-Scalar-Tensor theory [30 



(MST, where 



ao = d\nA{^o)/d(po, = 9a((^o)/5^o); (4) Scalar-Tensor- Vector theory [l| (STVG); (5) 
Modified Scalar- Tensor- Vector theory in this paper (MSTVG); (6) Einstein-aether theory 



48| (iE); (7) The tensor- vector-scalar theory provided by |49| (TeVeS, the results 7 and /? 
given by [5(j]). 

For different theories listed in Table [HI EEP is satisfied except MST, STVG and MSTVG. 
PPN formalism only works in validity of EEP so that comparisons with PPN parameters of 
all kinds of gravity theories, especially MST, STVG and MSTVG, is only phenomenological. 
We are going to argue that even though the PPN formalism is only valid under EEP, 7 
should not depart from 1 too much in STVG. Generally, the deflection of light can be used 
to test the parameter 7. Equations of motion for photons in IPN contain only two metric 
coefficients: N and V. For MSTVG, equations of motion for photons in IPN still contain 
the contribution of vector field besides the two metric coefficients. When considering only 
one point mass, equations of motion for photons in IPN for MSTVG is as follows 



GM ■ 
— 



GM , 
a — =-r 



7- 



GM 



dx 



dt 



/ + 2(l + 7)- 





dx* 







= g^,v- 



dt dt 

Let us assume the Newtonian order solution of these equation is 



(75) 
(76) 



x^ = cfCXt - to), 



(77) 
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where = 1 and light in the Newtonian order travels in a straight line at constant speed 
c. Furthermore, we assume the first order solution of these equation is 



x' = cn\t - to) + xi, (78) 



where is the initial emitting direction of a light signal. Then, substituting Eq. flTHl) into 
Eqs. dTSD-dZSl), we obtain 



d^Xp _ gM 



+2(l + 7)^n(n-r), (79) 
n-^ = -e(l + 7)^. (80) 



Followed the approach of Ref. 



24| . by using Eqs. fl79l)-fl80l). the light deflection angle up to 



IPN approximation for MSTVG is 

A</. = 26^^(1 + 7 + «e-'^/'), (81) 

where d represents the coordinate radius at the point of closest approach of the ray. Then 
the deviation from GR is 

5, ^ - = _ 1 + ^^-^lA ^ (82) 



which comes from two parts: PPN parameter 7 and the vector field. If 7 = 1/3, the case of 
Eq. fl8ip and fl82l) will return to STVG [l|. Here, one would claim that the part from the 
vector field contributes to remaining 2/3 which could make up 7 7^ 1 in STVG. However, this 
is not the case. When considering the Lense-Thirring drag on the orbit of an earth-satelite, 
the longitude of the ascending node f2 and the argument of pericenter lo vary with time by 



using the averaging method 



d^ ^ {l+l)ge^Jearth . ^ 

dt a3(l-e2)3/2 ' ^^"^^ 

du 2(1 + 7)^6^ Jearth . foA\ 

^ = - a3(l-e^)3/2 ''''^ 

where Jearth is the angular momentum of the earth and i is the inlination. Eqs. ( l83l) and 
(18^ return to GR when 7 = 1. In Lense-Thirring effect, it worthy of note that dQ/dt, 
namely Eq. (183|) . only depends on the force which is perpendicular to the orbital plane. In 
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STVG, however, the Yukawa force can not affect dQ/dt because it is a radial force in the 
orbital plane. For this reason, dfl/dt can limit 7 very closely even though STVG violate the 
EEP. This is the reason why we have to modified STVG. 

Now, we mainly focus on the three gravity theories that include a vector field. They are 
respectively IE TeVeS 49|] and STVG in this paper. iE is a tensor- vector theory in 
which the vector field is massless, unitary and time like. This theory investigates preferred 
frames and violation of Lorentz invariance. There are five gravitational and aether wave 
modes for iE-theory 5l(|. Recently, Xie and Huang 52] presented a 2PN approximation of 
JE-theoTj and found that the linearized waves with the spin-0 and spin-1 modes in JE will 
propagate with infinite velocities if the 2PN light deflection angle in iE-theory is identical 



53 



54 



55l | has gained recognition 



with that of GR. MOdified Newtonian Dynamics (MOND) 
as a successful scheme for explaining galaxy rotation curves without invoking dark matter. 
TeVeS is a scenario of relativistic MOND, which includes a unit massless timelike vector and 
a scalar field. In TeVeS, modified Newtonian acceleration is from two positive dimensionless 
parameters K and k but not scalar or vector field. TeVeS passes the usual solar system tests 
and provides a specific formalism for constructing cosmological models. For STVG, there 
are three scalar fields, one vector field of rest mass m besides the metric field. Both STVG 
and TeVeS try to explain galaxy rotation curves without dark matter, but the modified 
Newtonian acceleration in STVG is from the vector field. 

For MSTVG in IPN, other eight PPN parameters except 7 and (3 are all zero. Measure- 
ment of the defiection of light passin g Ju piter by Very Long Baseline Interferometer (VLBI) 
gives 7 — 1 = (—1.7 ± 4.5) x 10^^ 56|. The most precise measurement of 7 comes from 
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the Cassini experiment by Doppler tracking, which gives 7 — 1 = (2.1 ± 2.3) x 10^^ 
By comparing the masses of 15 elliptical lensing galaxies from Sloan Lens Advanced Cam- 
era for Surveys (SLAGS) on kiloparsec scales (about lO^^m), Bolton et al. 57| constrain 
7 = 0.98 ± 0.07 in Icr confidence level (GL). For (3, although its accuracy level is lower than 
that of 7, some experiments give the limits of jS. For example, /3 — 1 = 3 x 10^'^ comes 



from t 
effect 



le perihelion shift of Mercury 58|] and (3 



2.3 X 10 comes from the Nordtvedt 



59 



60l |. Tests of 7 and (3 are basically going on the solar system scale. However, 



according to present-day experiment data, the values of these two parameters (7 and f3) are 
independent of the scale of the testing systems. In order to obtain 7 = 1 and (3 = 1, the 
parameter ^0 in MSTVG must go to infinity with 61 growing slower than ^q. 
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TABLE III: Yukawa parameters a and A in macroscopical range. 





\ol\ 


A 


Ref. 


SDSS 


0.35 ± 0.9 


10Mpc/h< A <100Mpc/h 


[63.1 


2dFGRS 


0.025 ± 1.7 


10Mpc/h< A <100Mpc/h 


[631 


Pioneer anomaly 10/11, 








Galileo, and Ulysses data 


10-3 


200AU 


[691 


Long range limit 


> 10-4 


> 70AU 


[701 


Pioneer anomaly fSTVG) 


10-3 


47 ± lAU 


[111 


Constraint from Sun (STVG) 


10-8 


lO^^m 


[151 


Constraint from Earth (STVG) 


10-13 


10i3m 


[151 


Planetary data(EPM2004) 


IQ-i^ ~ 10^13 


lAU 


[651 


Constraint on solar system 


(0.3 ± 2.7) X 10^11 


0.2 ± 0.4AU 


[4Q1 


planetary motions 


< 10-9 


< 0.18AU 


[661 


T T T3 


^ in-io 

> iU 


lu ni<, A lu ni 




planetary data 


> 10-9 


> lOV 




Constraint on intermediate-range gravity 


« 10-8 


1.2 X 10'^m< A < 3.8 X 108m 


m 


LAGEOS-lunar 


10-8 < \a\ < 10-5 


lO^nK A < lOV 


[36] 


Earth-LAGEOS 


10-*^ < \a\ < 10-3 


lO-^nK A < lO^m 


[36] 



B. Review and discussion on Yukawa parameters 

rn 

Since Fischbach et aL [61j] suggest a possible gravity-like "fifth" fundamental force in 
macroscopic scale, it evokes some interest in many theories which intend to unify gravity 
with other known forces. The presence of the fifth force could be detected by searching 
for apparent deviations from Newtonian gravity. For instance, the fifth force would arise 
from the exchange of a new ultra-light boson which coupled to ordinary matter with a 
strength comparable to gravity. Typically, through added a hypothetical Yukawa force to 
the Newtonian potential, this modified potential per mass takes the form: 

[/(^) = _^(^l + «e-^/A^, (85) 

where a represents the strength of the Yukawa coupling, and A represents its length scale. 
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In Table llllt we list some resu 
10"^m. Fischbach and Talmadge 



ts about a and A in macroscopical range from lOOMpc to 



36j consider the model of Eq. fl85|) in astronomical tests 



which provide tight constraints on Yukawa parameters [a and A). Typically, these results 
are based on testing G{r)MQ values deduced for different planets through the following 
equation 



G{r)M 



where 



G(r) 



l + 1 + ^ le 



-r/A 



(86) 



(87) 



Planetary data gives |a| > 10 ^ and A > lO^m (see Table UTTjl . The constraints at larger 
ranges from laboratory, geo physical, and astronomical data are essentially unchanged (for 



detail, see Figure 2.13 of [36|]). As to the constraint from LLR, Adelberger et al.[62| updated 
the result to include recent LLR data and gave |a| > 10^^° and 10^ < A < lO^^m. 

By using two large-scale structure surveys: the Sloan Digital Sky Survey (SDSS) and 
;wo-degree Field Galaxy Redshift Surveys (2dFGRS) on megaparsec scales, Sealfon et 



the 
al. 



63| considered two models about Poission equation which deviated from gravitational 



inverse- square law to constraint a with marginalized over A from lOMpc/h to lOOMpc/h. 
Where h is the value of Hubble constant in units of lOOkm/s/Mpc. The potential of the 
first model is 



d'^r 



1 + a 1 



-|r-r'|/A 



(88) 



After integration of Eq. (IHHI) . it has the same result as the Eq. (l85l) . On this large-scale 
structure scales, a = 0.025±1.7 for 2dFGRS and a = — 0.35±0.9 for SDSS at 68% confidence 
level (CL) through fitting the power spectrum measurements from SDSS and 2dFGRS. 

Li and Zhao [64] consider Eqs. fl5Bl) - fl57l) and constrain the Yukawa parameters. Using 
the earth-satellite measurement of earth gravity, the lunar orbiter measurement of lunar 
gravity, and lunar laser ranging measurement ( 6J| ) , the result from constraint on the two 
Yukawa parameters are a = 10"^ — 5 x 10~^ and A = 1.2 x 10''m— 3.8 x lO^m. 

Moffat and Toth [L5\ express the acceleration law in STVG as 



where 



G. 



1 + a - a ( 1 + - 1 e""/^ 



(89) 



(90) 
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With observations or experiments performed within the solar system or in Earthbound 



laboratories, the authors estimate A = 5 x 10 m, |a| = 3 x 10 for the Sun, and A = 
8.7 X lO^^m, |a| = 9 X 10"^^ for the Earth. 

By analysis of EPM2004 ephemerides, lorio Gsl constrains on the strength and rang of a 



and gives |a| 



10 



-12 



10 



-13 



66( 1 constrains on the range and the 



Yukawa-like fifth force through considering potential Eq^ 
A ~ lAU. By using the same potential Eq. (!85l) . lorio 
strength of a Yukawa-like fifth force with planetary perihelia by the EPM2004 ephemerides: 
A < 0.18AU and |a| < 10"^ Through considering Eqs. ([89D-([90D, lorio |40| obtains A = 
0.2 ± 0.4AU and |a| = (0.3 ± 2.7) x 10^^^ by constraint on perihelion shift of Mercury. 

With Pioneer 10/11 launched, their radiometric tracking data have consistently indicated 
the presence of a small anomaly which is = (8.74 ± 1.33) x 10^^'^m/s^, directed toward 
the Sun This apparent anomalous constant acceleration has been similarly shown in 
Galileo and Ulysses data 68|]. In order to explain this anomaly, some of recent efforts are 
looking for new physics. Based on this idea, many modified Newton inverse square laws were 
provided. John et al. 69[ ruled out a lot of potential causes and considered a gravitational 
potential with introducing an additional Yukawa force but this modified potential is a little 
different from Eq. fl85|) and has the following form 

QM 



V{r) 



ae 



-r/X 



(91) 



Reynaud and Jaekel 



(1 + a)r 

ast term of Eq. fl9Tl) as the Pioneer 10/11, Galileo and Ulysses 



Through identifying the 
acceleration, John et al. 69 1 obtain a = — 1 x 10^^ for A = 200AU. 



70| also discussed the relation between long range tests of the 



Newton law and the anomaly recorded on Pioneer 10/11 probes through considering the 
potential of Eq. (jHSj). With a power expansion of Eq. fl55]) in terms of r/X, Reynaud and 
Jaekel [t^ obtain a constant anomalous acceleration on the range of distances probed by 
Pioneer 10/11 only if A > lOAU (a = -AVA^ and A = 6300AU). It yields a < -10"^ 
or |a| > 10~^. In STVG, Brownstein and Moffat 7l| considered the anomaly recorded 
on Pioneer 10/11 probes through Eqs. ([89])-([90D and obtained a = (1.00 ± 0.02) x 10^^ 
A = 47 ± lAU. For these results, it can be seen that a long-range Yukawa deviation from 
the Newton potential can be treated as a constant acceleration. 

From above analysis, the Yukawa parameters depend on the scale of testing system. In 
galactic scale, Yukawa parameters would play a big role to explain flat rotation curve Q]. 
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TABLE IV: Parameters in binary pulsars. 



Pulsars e M (M©) P (day) < Co > (°yr~i) Ref. 

PSR B1913+16 0.6171338(4) 2.8281(2) 0.322997462727(5) 4.226595(5) [76] 
PSR B1534+12 0.2736767(1) 2.679(2) 0.420737299153(4) 1.755805(3) [72] 
PSR J0737-3039 0.0877775(9) 2.58708(16) 0.10225156248(5) 16.89947(68) [75] 
PSR B2127+11C 0.681395(2) 2.71279(13) 0.33528204828(5) 4.4644(1) [74] 

However, in the outer solar system, Yukawa parameters could provide to a certain extent 
the Pioneer anomalous constant acceleration. Then, in the inner solar system, the strength 
coupled with gravitation about the Yukawa force, a, is very small and the range. A, is large. 
These results lead to the fifth force negligible in inner solar system. If future experiments 
confirm nonentity of dark matter in the universe, it will provide a possible existence of a 
fifth force. Until then, the fifth force must work in galactic scale even cosmological scale. 
From this point of view, constraint on the Yukawa parameters [a and A) in MSTVG has 
somewhat practical significance. 

In view of these discussion about the four parameters in MSTVG: PPN parameters (7, 
j3) and the Yukawa parameters (a. A), the two PPN parameters are independence of the 
system scale in the current experiments. However, the case for the two Yukawa parameters 
are not. In the next section, the orbital data of pulsar binaries are used to fit the Yukawa 
parameters, in which the PPN parameters will be fixed as 7 = /3 = 1. 



IV. CONSTRAINTS ON THE YUKAWA PARAMETERS BY DOUBLE NEU- 
TRON STAR BINARIES IN MSTVG 

Double neutron star binaries provide us more impressive tests of general relativity than 
other systems. For example the fraction which will merge due to gravitational wave emission 
is larger. Besides, a neutron star is rather compact and its companion hardly effect its 
shape. These systems are highly valuable for measuring the effects of gravity and testing 
gravitational theories. Four samples used in this paper are PSR B1913+16, PSR B1534+12, 
PSR J0737-3039 and PSR B2127+11C which are listed in Table lYl The numbers inside a 
pair of parentheses are the la error of its corresponding quantity. 



PSR B1913+16 is discovered in 1974 by using the Arecibo 305m antenna ISj]. The orbit 
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5 10 15 20 25 30 

?./lonm) 



FIG. 1: The plane of \-a constrained by lo of four double neutron star binaries, where the real 
line (green) denotes PSR B1913+16, the thinnest dashed line (red) denotes PSR B1534+12, the 
thickest dashed line (blue) denotes PSR J0737-3039, and the dotted line (magenta) denotes PSR 
B2127+11C. It can be seen that the four curves are nearly cross at one point. 




\ PSR B21 27+1 IC 

""'"--^^ PSRB1534+12 
4 5 6 7 8 9 

A./10Hm) 

FIG. 2: The enlarged diagram of the region around the cross point in Figure [H The legend is the 
same as in Figure [TJ 
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1.5- 

3.6 3.8 4 4.2 4.4 4.6 4.8 

A./io'(m) 

FIG. 3: Each curve in Figl2]is separated into two curves to show the bounds due to the observational 
la error. Please read the context for the details. Figure legend is the same as in Figured! 
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A./10'(m) 



FIG. 4: Confidence level contours of 68.3% (red), 95.4% (green) and 99.73% (blue) in the plane 
of X-a constrained by four double neutron star binaries. The plus sign denotes the 68.3% values 
for the Yukawa parameters of the MSTVG scenario which are: A = (3.97 it 0.01) x 10®m and 
a = (2.40 ± 0.02) x 10~^, the minimum value of is 21.47. 
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has evolved since the binary system was initially discovered, in precise agreement with the 
loss of energy due to gravitational wave emission predicted by Einstein's General Theory 
of Relativity. PSR B2127+11C is located in the globular cluster Messier 15. This system 
appears to be a clone of PSR B1913+16 (see Table Mi)- For PSR B1534+12, its pulses 
are significantly stronger and narrower than those of PSR B1913+16 Q- PSR J0737-3039 
was discovered during a pulsar search carried out using a multibeam receiver with 64m 



radio telescope in 2003 



73| . And we can see this system has shorter orbital period, smaller 



eccentricity and larger periastron advance. When we fix 7 = 1 and /3 = 1, each u in Table 
HVlwill confine the values of a and A in a curve through Eg. ( 1731) (see Fig. [Hand Fig. [2]). 

In Fig. [H the trajectories of four double neutron star binaries in the plane of the Yukawa 
parameters have been shown. From Fig. [H we can see that the evolution of a — A by PSR 
B1913+16 and PSR B2127+11C is very similar. For PSR B1534+12, the evolution of a - A 
is mostly flat. Besides, the evolution of a — A for PSR J0737-3039 is very steep. In Fig. 
[H we can also see that the four trajectories almost meet together at ~ 4. Fig. [2] is 

the enlarged drawing of Fig. [T] at a smaller scale. On the other hand, Damour & Esposito- 



Farese 



301] constraint two 2PN parameters of MST by using four different binary pulsars 



data in la confidence level. Based on the same method of [30], we constrain the Yukawa 
parameters of MSTVG in la confidence level. Firstly, we plot an la constraint imposed by 
double neutron star binaries. Each set of binary data leads to a reduced x^'- xlinaryi^^ ^) — 

[y^theory ^observation )^ / (^observations equivalent to the la constraint -aobservation < ^theory - 

observation < observation- The bouuds by four double ueutrou star binaries allowed regions 
of the A-a plane are displayed in Fig. [3l Clearly four binaries data favor only a small region 
of the Yukawa parameters. To combine the constraints on a and A coming from different 
double neutron star binaries experiments, we have added their individual as if they were 
part of a total experiment with uncorrelated Gaussian errors like the analysis method of 30| : 
xLaz(«, ^) = X?9i3+i6(«> ^) + X?534+i2(«> ^) + Xo737-3039 ("> ^) + xii27+iic;("> Therefore, 
we plot the contour level AxLa;(". ^) = 2-3, AxLa«("' ^) = 6-17, and AxLa«("' ^) = H-S^ 
where AxLa/(">^) = xLai("'^) - (xLai(«> ^))mm, defines respectively for two degrees of 
freedom the 68.3%, 95.4% and 99.73% confidence levels represented in Fig. |H The la fit 
values for the Yukawa parameters are A = (3.97 ±0.01) x lO^m and a = (2.40 ±0.02) x 10~^ 
with Xmin — 21.47. With these results, we go back to discuss the light deflection for Eq. 
(1821) . If we consider that a light ray which passes the Sun at the solar radius, it yields 
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6^ = 3.33 X 10~^. This means the deviation of MSTVG from GR for defection of hght may 
should be tested in the future accuracy level. 

For the same model of the Newtonian potential modified by the fifth force (Eq. fl85|)). we 



obtain the same results as in LLR 
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62| , constraint on intermediate-range gravity [64] and 



LAGEOS-lunar [36] listed in Table [XTTl It tells us that the limits of binary pulsar systems 
on the Yukawa parameters for MSTVG are basically consistent with the solar system. From 
the other view, it tells us that |27 — /3 — 1| almost equals to 0. When considering three 
parameters {a, A and \2'y — (3 — 1\) based on Eq. (!73|) in la level, we do not obtain any 
reasonable results and need more precision binary pulsars data. 



V. CONCLUSIONS AND PROSPECTS 



With relaxed the two assumptions taken in STVG, which actually hold when the scalar 
field G is a constant field, it is shown that the post-Newtonian parameter 7 7^ 1 by using 
Chandrasekhar's approach and a modified scheme of scalar-tensor-vector gravity theory 
(MSTVG) is then proposed by introducing a coupling function of the scalar field G: 6{G). 
Started with the action in MSTVG, the equations of motion of MSTVG in the first post- 
Newtonian order (IPN) for general matter without specific equation of state and N point 
masses are obtained. The secular periastron shift for binary pulsar in IPN is derived. 
From the results of MSTVG in IPN, there are four parameters: two PPN parameters 7 
and (3 and two Yukawa parameters a and A. After pointed out their independence of 
system scale for 7 and /5, discussion about the dependence of system scale for the Yukawa 
parameters is touched. Furthermore, with the precondition of |27 — /? — 1| = 0, a and A in 
MSTVG with applied 4 double neutron star binaries data (PSR B1913+16, PSR B1534+12, 
PSR J0737-3039, PSR B2127+11C) are constrained: A = (3.97 ± 0.01) x lO^m and a = 
(2.40 ±0.02) X 10"*^. It has been shown that the limits of binary pulsars systems on MSTVG 
are basically consistent with the results from the solar system such as the earth-satellite 
measurement of earth gravity, the lunar orbiter measurement of lunar gravity, and lunar laser 
ranging measurement to constrain the fifth force. For future applications in binary pulsars 
systems, it can help us to distinguish between different gravity theories and MSTVG. If 
MSTVG is proved to correct, besides success in solar system, astrophysical and cosmo logical 
scales without dark matter, there must exist a wave of vector field which does not equal to 
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the light speed. Besides, there also exist violations of the Einstein equivalence principle in 
large scale in future experiments. 

In Appendix B, we calculate the standard PPN parameters of the MSTVG, ignoring the 
vector field due to the coupling between it and the matter fields. When the vector field is 
included, some new super potentials might be introduced that would cause the appearance 
of new parameters and the numerical values of some PPN parameters might be affected. 
This is a subject of future research. 
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APPENDIX A: CALCULATION OF 7 

In this appendix, we will give a detailed calculation of 7 in MSTVG by using Chan- 
drasekhar's approach. The contravariant components of the metric tensor are 



We can now evaluate the Christoffel symbols with the aid of the metric coefficients given 




(Al) 



(A2) 
(A3) 



Accordingly, 




(A4) 
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in Eqs. 



and flXB- flM]) : 



^ 00 



,1 
2" 



^ 00 



Then, the general expression for the Ricci tensor is 



Rij 



1 



1 



-^\N^kk 



\N^kN,k - ^H^tt + ^N^ikHik + ^N^kHki,i 



— A'^ i H k — — L kk ~\~ Lk kt 



Turning to the components of the energy-momentum tensor, we find 



^00 


= cr - e 




= -ec^i, 


T- 


2 



and 



T 



kk 



The (0,0) component for from Eq. (fT2|) : 



- -UN 




2 





1 {2)- 

A-nG^a + -DC 



The component for Hij from Eq. (fT2|) : 



— -Hij^kk + -Hik.ki + -Hjk.ki — -H, 



ik,kj 



hjk.ki 



(2) 



(A5) 
(A6) 
(A7) 
(A8) 
(A9) 
(AlO) 

(All) 
(A12) 



(2)1 



(A13) 

(A14) 
(A15) 
(A16) 

(A17) 
(A18) 

(A19) 



29 



Using the gauge equation (l32l) . we obtain 



2 Hij^kk 



1 {2)' 



(2) 



For G from Eq. (|20l) 



(2)- 

9o + 3)nG' 



Based on Eqs. (lAlSP and (1A21|) . we obtain 

^0 



^0 + 3 



SttGo. 



Based on Eqs. f!A2np and f!A21l) . we obtain 



If we define if,. 



5ijV , then we obtain 



and 



UV 



1 



^0 + 2 
^0 + 3 

^0 + 2 
^0 + 3 

V ^0 + 2 



SvrGo, 



(A20) 



(A21) 



(A22) 



(A23) 



(A24) 



(A25) 



A^ do + A 

When 9q = —1, it returns to the case of STVG and gives 7 = 1/3. It is also exphcit that 

(2) 

7 = 1 when the scalar field G is a constant field. In this special case the terms with G 
in Eqs. f lATSl) and flA20l) disappear and 7 = 1 holds for both STVG and MSTVG. One 

(2) 

might wonder at Eq. (1A21I) . which does not allow G to be zero. But Eq. ( 1A21I) comes from 
Eq. (I20!) . the field equation for the scaler field G: Eq. (!20l) does not exist when G is not a 
variable field. 

The next step is to identify the parameter 7 defined in Eq. (1A25P as the corresponding 
PPN parameter. This is done in Appendix B. 



APPENDIX B: DERIVATION OF THE PPN PARAMETERS FOR MSTVG 

In this appendix, we mainly derive the PPN parameters of MSTVG. For this purpose, 
we must transfer our coordinate into the standard PPN coordinate and then derive PPN 
parameters of MSTVG in comparison with PPN formalism {24]. But, it is shown that 
MSTVG violates EEP based on Eq. ([8]) and can not be brought into the standard PPN 
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formalism. To compare with the PPN metric, however, we provisionally ignore the vector 
field of MSTVG, which causes the violation of EEP. 

We consider the material composing the various bodies of the system to behave like an 
ideal fiuid. Following T^], the energy- momentum tensor can be written in the following 
form in the ideal fiuid case 



(Bl) 



where p denotes the rest-mass density, p is the pressure, 11 is the specific internal energy, 
and = dx'^/cdr is the dimensionless 4- velocity and we obtain in IPN 

(B2) 

(B3) 
(B4) 

where is the coordinate velocity of the corresponding material element. 
From Eqs. ([I1D-([I5]) and ([12]) -([£4]), we obtain 



rpOO 


= P 


l + \{U + v^ + N) 






c 




rpij 


= ^{pv'v^ + p6ij) + 0{c 





a = p 



1 + -(n + 2v^ + N) 



+ 34 + 0{c 



= pv' + Oic-^), 



(B5) 
(B6) 
(B7) 



Then, we rewrite Eqs. (l35l) and ( l42l) with abandoning the vector field 



N + e'L 



-SnGp 

V ^0 + 4;^ p^ 

+N 



+ 2t;2 + n + 3- 



3^iV-2(l 



2(^o + 3)(^o + 4)^ 



N-N 



2V ^2(^^o + 3)(^^o + 4)2 



AN' 



tt f, 



m 
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Then, we obtain 
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With the above definition of the gravitational potentials, it yields solution of the metric for 
MSTVG without the vector field in the following forms 



N 
L 



2U, 

-2(1 + 
+2 



(B25) 



01 



2(^0 + 3)(^o + 4)2 



f/2 + (2 



00 + 2 
^0 + 4 



+ 2)$i 



'00 + 2 / 01 
.^^;;+4 ~ + 2(^0 + 3)(^o + 4)2 



+ 1 



$2 + 2<l'3 



+6 



^0 + 4 



-$4 



X,tt, 



u 


= -2(^1 + 




00 + 2 
= '00 + ^' 



00 + 2 

00 + 4: 



(B26) 
(B27) 
(B28) 



In order to obtain PPN parameters, we must transfer our coordinate system into the 
standard PPN one. And then we could compare the metric of MSTVG with the PPN metric 
in the standard PPN coordinate system and finally derive the PPN parameters of MSTVG. 
When the above is compared with the standard PPN metric (see Eqs. flB43P - flB45l) . the only 
superpotential which does not appear in the PPN metric is x,tt in Eq. (]B26p . Therefore, it 
is necessary to transform the coordinates to gauge off this term. Based on the Eq. (4.38) 
in Ref 2^, an infinitesimal coordinate transformation is introduced between the standard 
PPN coordinate system and ours: 



where 

^0 = AiX,o, = A2X,i- (B30) 
The relation between the metrics before and after the gauge transformation, g^i, and g^u 



respectively, are shown by Eq. (4.46) in Ref 

- _ 2 \ 

9ij — 9ij ~ f ^sX.jj) 

9oi = 9oi - e^(Ai + A2)x,ti, 

goo - e''2\ix,tt - e'2A2(t/' + $^ - $2) 



(B31) 
(B32) 
(B33) 



Due to the spatial part of the PPN metric and our metric are all diagonal and isotropic, we 
thus choose A2 = through Eqs. flB3Tl) and substitute Eq. flB28l) into Eqs. flB3B : 



9ij 



l + e'2^-^U\5., 



00 + 4 



(B34) 
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For Eqs. flB32l) and flB33|) . we have 



-e"2 1 



^0 + 4 



(B35) 



goo = goo - e^2Aix,« 
= -1 + e^2U + 



-2 1 + 



+2 



^0 + 2 
3t^^^ - 21 1 



^0 + 4 
-(2Ai + 



2(^0 + 3)(^^o + 4)2 



, ^0 + 2 
'^^0 + 4 



+ 2 $1 



$2 + 2<l>3 + 6^^<l>4 

^0 + 4 



(B36) 



by using Eqs. (1B25|) . flB26|) and (1B27|) . It is noted that there is no existence of superpotential 
X,tt in the standard PPN framework, we then obtain the following value through Eq. flB36P 



Ai 



1 

2 



(B37) 



When we substitute A2 = and Ai = — | into Eq. flB29p . the transformation between our 
reference system and the standard PPN reference system is shown as 



tpN = t + e^-x,t + 0{5), 



X 



PN 



(B38) 
(B39) 



Through this transformation, our metric in the PPN coordinate system become 



gij 



1 + e^2 



Qo 



o 1 / 6^0 + 2 \ o 1 
^00 = -l + e22f/ + 6^1-2^1 + 



+2 



^0 + 2 
"i^- 7-2 1 + 



2(^0 + 3)(0o + 4)2 



+ 1 



(B40) 
(B41) 

^^+(2^ + 2)$, 
V ^0 + 4 ) 

<l>2 + 2<l>3 + 6^^<l>4j, (B42) 

[70 + 4 I 



■^0 + 4 "V^ ' 2(^0 + 3)(^o + 4)2^ 
by using Eqs. (lB34D - (lB36|) . 

On the other hand, with the above definition of the gravitational potentials, the standard 
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PPN metric 2J] reads 



9ij 
900 



(1 + e^2^U)5, 

^3 



i(47 + 3 + ai - «2 + Ci - 20Vi - ^(1 + ^2 - Ci + ^^W, 



(B43) 
(B44) 



-1 + e^2U + - 2(3U^ - 2^$^ + (27 + 2 + + Ci - 20$i 
+2(37 - 2/3 + 1 + C2 + 0*2 + 2(1 + C3)$3 + 2(37 + 3C4 - 2e)$4 
-(Ci - 2021 



(B45) 

Thus, the PPN parameters of MSTVG without the vector field have the following forms by 
comparison between Eqs. flB40|) - (lB42l) and (iB43|l - flB45l) 

O0 + 2 



7 



P 



^0 + 4' 
1 



01 



2(^^o + 3)(^o + 4)2' 

C = "1 = "2 = "3 = Cl = C2 = Cs = C4 



(B46) 

(B47) 
(B48) 
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